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It is argued that, if a regular Hamiltonian is perturbed by a term that produces chaos, the onset 
of chaos is shifted towards larger values of the perturbation parameter if the unperturbed spectrum 
is degenerate and the lifting of the degeneracy is of second order in this parameter. The argument 
is based on the behaviour of the exceptional points of the full problem. 
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I. INTRODUCTION 

The study of single particle motion in deformed mean fields has attracted much attention recently because of its 
relevance to nuclear physics and for the description of metallic clusters ||,^. For the case of harmonic oscillator 
potentials, deformations which go beyond a quadrupole deformation are nonintegrable and show chaotic behaviour in 
the classical and quantum mechanical treatment. However, closer scrutiny has revealed that the addition of an octupole 
term to a prolate quadrupole potential produces the typical signatures of chaos only for fairly large octupole strength, 
in fact, the problem appears to be close to integrability. Moreover, the quantum mechanical treatment produces new 
shell structures for finite octupole strength even though the classical problem is nonlinear |^,^. In contrast, addition 
of an octupole term to an oblate quadrupole deformed potential yields chaos with a positive Lyapunov exponent, and 
the corresponding quantum spectrum has the typical statistical properties ascribed to quantum chaos. 

In this paper we address the question: what are the intrinsic properties of the quantum mechanical operators which 
give rise to the different behaviour described above. This question is of interest for a possible characterisation of what 
is called quantum chaos. The example mentioned in the previous paragraph renders a good case for such studies as 
it refers to the orthodox situation where a classically chaotic system is treated quantum mechanically. Our aim is 
to unravel the universal operator properties that produce the typical patterns ascribed to quantum chaos without 
reference to an underlying classical system. There is progress in previous work towards this aim as it has been 
recognised that, for a problem of the form + XH^, it is precisely the high density of the exceptional points [^|j, 
i.e. the singularities of the spectrum En{X), which bring about the statistical properties of the spectrum associated 
with quantum chaos . The new aspect in the present paper is the effect of degeneracies at A = 0 upon the behaviour 
for A > 0. We find that if the lifting of degeneracies is a second order effect in A, then the onset of chaos is suppressed 
and will occur only for larger values of A. It is this pattern which explains the difference between the prolate and 
oblate case as described above. The major point of the present paper is the universal validity of this finding, and the 
behaviour of the exceptional points around A = 0 provides for the explanation. 

In the following section we introduce briefly the concept of exceptional points and present the argument for the 
statement made above. Section three contains examples for illustration and section four refers to the particular 
physical example introduced above. A summary and discussion is given in the last section. 


II. EFFECT OF DEGENERACIES ON EXCEPTIONAL POINTS 

To make the paper selfcontained we briefly recapitulate the significance of exceptional points and their connection 
to avoided level crossings. This will facilitate the discussion about the effect of degeneracies on exceptional points 
and hence the global structure of the spectrum. 

There is essentially a one to one relationship between avoided level crossings and exceptional points [||. If we 
have a quantum mechanical problem of the form + XH^ with and given as N x N matrices, then the 
spctrum En(X) ,n = 1,... ,N is determined by one analytic function evaluated on N Riemann sheets 1^. If there 
are no degeneracies, the N sheets are connected by N{N — 1) branch points, the exceptional points. They occur in 
complex conjugate pairs and are the points in the complex A-plane where any two pairs of energies coalesce. If this 
happens sufficiently close to the real A-axis, an avoided level crossing occurs since the pair of energies coalescing in the 
complex plane still assumes values near to each other on the neighbouring real A-axis. In principle, the positions of the 
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exceptional points are determined by the resultant of the secular equation for the spectrum, i.e. by the simultaneous 
solution of the two ploynomial equations 

- \H'^) = 0 

^det(i;-7J“-Ai?i) =0 (1) 

which leads by elimination of E to the resultant which is a polynomial of order N{N — 1) in A. For large values of N, 
an explicit calculation of the exceptional points is prohibitive. However, it is possible to determine a distribution for 
the real parts of the exceptional points from the knowledge of and iJ^ alone , and in this way the regions of real 
A-values with large scale avoided level crossings can be found. These are the regions where the spectrum shows the 
properties ascribed to quantum chaos. This method has been tested and applied to the particular physical situation 
of the hydrogen atom in a strong magnetic field [Q. 

The effect of removing degeneracies by small random perturbation was discussed previously In the present 
paper we address the particular problem of lifting degeneracies by the term Ai/^, i.e. we assume that El'^ has certain 
systematic degeneracies (like for example the harmonic oscillator). Guided by the particular physical problem as 
discussed in the introduction and resumed in section four, we consider and compare two situations which turn out to 
be significantly different in the context of quantum chaotic behaviour. The two situations are characterised by the 
presence or absence of a linear term in the expansion of En{X) in powers of A around A = 0. 

It is obvious that the spectrum has a distinctly different behaviour in the vicinity of A = 0 for the two cases 
considered. A non-vanishing linear term exhibits the typical behaviour of a lifting of degeneracies, when a perturbation 
is switched on, in that the levels fan out of the degenerate level when A is turned on. With a zero linear term the 
levels stay closely together and separate only when the quadratic term becomes significant. From this picture one 
intuitively expects that, if the additional term \E[^ produces chaos at all, the onset of chaos is delayed for increasing 
A in the latter case when comparing with the former case. In the following we confirm this expectation using the 
exceptional points for the argument, while the next section presents some illustrative examples. 

Each m-fold degeneracy reduces the total number of exceptional points by the number m{m— 1) if the linear terms 
are present in the expansion of the energy levels. If the leading perturbing terms are quadratic, there is a further 
reduction by the same number m{m — 1). This follows from the resultant which starts with as the lowest 

order in the former case and with ;^ 2 m(m-i) latter case. This behaviour is valid for each m-fold degeneracy, 

in other words, if the levels have at A = 0 the degeneracies mi, m 2 ,..., the lowest order terms of the resultant start 
with i-jjg gf |.j^g degeneracy is linear and with ... foj. 

a quadratic lifting. A particular m-fold degeneracy can be viewed as a confluent situation where m(m — 1) complex 
conjugate branch points have merged into the point A = 0 thus cancelling the singularities altogether; likewise, in 
the quadratic case, 2m(m — 1) singularities have merged. The essential point of the argument lies in the comparison 
between the two cases: switching on linear terms leads to the emergence of — 1) additional exceptional 

points. The distance from A = 0 of the additional exceptional points depends on the magnitude of the linear terms. 
For typical values they are rather closer to than remote from the origin. For large matrices the additional number of 
exceptional points can amount to many thousands. This causes avoided level crossings in the spectrum in a region 
where without the linear term the spectrum appears smooth. 


III. ILLUSTRATION OF SIMPLE EXAMPLES 

The realistic example which initiated this work produces more than twenty nine thousand exceptional points for the 
matrix size used in our work . We return to its treatment in the following section. Here we illustrate the behaviour 
of the exceptional points in a low dimensional example. The effect on the spectrum, in particular the different onset 
of quantum chaos for the two different cases, is subsequently demonstrated in a generic matrix model. 


We consider a nine dimensional model where E[^ has the 

diagonal 

entries (1, 

1,1,2,2,2,3,3,3). The three-fold degen- 

eracies are lifted by the term XEl^ with 
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For n = 0, has only zero entries in the degenerate subspaces of therefore the nine levels of the full problem 
have a vanishing linear term when expanded at A = 0. Otherwise the entries of have no significance, they have 
been generated randomly, any other choice (with similar orders of magnitude) yields the same qualitative result. Also, 
a diagonal form of in the degenerate subspaces can always be achieved by orthogonal transformations within the 
subspaces without changing the global spectrum. The specific form chosen is convenient as it leaves out additional 
irrelevant parameters and it is in line with the physical situation discussed in the following section. 

A nine dimensional model yields 36 exceptional points in each half plane, our particular case reduces this number 
to 18. In Fig.(la) we display most of the exceptional points in the upper A-plane. When the parameter /i is switched 
on, nine additional exceptional points emerge from A = 0 in each half plane. This is illustrated in Fig.(lb). The basic 
difference between the two figures lies in the additional singularities scattered around A = 0. Note that the other 
exceptional points change only moderately under variation of fi. The corresponding spectra are displayed in Figs.(2). 
The growing distance from the origin of the avoided crossings with increasing value ^ is clearly discernible in Figs.(2c) 
and (2d), where only part of the spectrum is shown for better illustration. These avoided crossings are due to the 
’new’ exceptional points. Note that in particular the levels which originate from the same degenerate energy are 
affected. When this occurs on a large scale, the nearest neighbour distribution of the spectrum will typically assume 
the Wigner surmise for rather small values of A while there is no resemblence to the Wigner curve when the linear 
terms of the perturbation vanish. This is demonstrated in the following study case. 

Essentially we repeat the model used above. We consider N > 100 for the full dimension and the dimensions of the 
degenerate subspaces can be two, three, four and so on. The diagonal matrix contains the entries (1,1,...,2,2,...) 
and the entries of are filled randomly with elements between —1 and +1. The off-diagonal elements of the diagonal 
blocks of the degenerate subspaces are set equal to zero while the corresponding diagonals are filled with the numbers 
/^(—’^deg/S + 1/2, —ndeg/2 -I- 1,..., ndeg/2 — 1/2) with n^eg being the dimension of the degenerate subspace. Because 
of the huge number of exceptional points, we no longer focus our attention on their positions. Instead, we compare 
nearest neighbour distributions for ^ = 0 and p. 7 ^ 0 for small values of A. 

In Figs. (3) typical results for the two different cases are presented. For the same value of A sample averages of 
ten samples of nearest neighbour distributions are illustrated for ^ = 0 and ^ = 1. There is of course a continuous 
transition from Fig.(3a) to Fig.(3b), and only if the order of magnitude of the diagonal elements of has reached 
that of the other matrix elements, that is for p, ~ 1, the Wigner distribution has fully developed. We can understand, 
in terms of the exceptional points, why for 0 < ^ <C 1 the situation closely resembles that of ^ = 0 even though a 
large number of exceptional points has already emerged from the origin of the A-plane. When the exceptional points 
are still very close to the origin, that is very close to each other, there is a cancellation with regard to the effect 
on the spectrum. In fact, two square root branch points connecting the same Riemann sheets, are barely noticeable 
from a distance; the function + looks just like for \z\ ^ e. Only when the branch points have moved out 

sufficiently far does the effect on the spectrum become significant as clearly seen in Figs. (3). The results presented 
refer to rideg = 8 and N = 304, which yields 2128 exceptional points having emerged from the origin. When rideg 
is increased, the window of A-values, where this transition is strongly pronounced, becomes wider. Conversely, for 
smaller values of rideg the effect is clearly seen only in a smaller range of A-values. This follows from the higher 
density of exceptional points associated with larger values of rideg- Also, the A-values beyond which a plain Wigner 
distribution occurs even when fi = 0 moves closer towards the origin with increasing rideg- These two observations can 
qualitatively be understood from perturbative arguments in that we expect the onset of a Wigner distribution beyond 
the intersection point of two straight lines or parabolae from neighbouring levels for fj, ^ 0 or fj, = 0, respectively; the 
curves intersect at Aintsct which is proportional to l/rideg in the former and l/^rideg in the latter case. This argument 
indicates where perturbation breaks down; we recall that the present paper deals with a nonperturbative situation as 
it is just the exceptional points which break down perturbation. We note that in a typical physical situation, like the 
one discussed in the next section, one is usually faced with steadily increasing values of rideg with increasing energy. 


IV. A PHYSICAL EXAMPLE 


Here we present the physical example which actually initiated this work. Phenomenological mean fields which 
contain quadrupole and octupole deformations have been investigated classically and quantum mechanically [QH] 
as this is of interest in nuclear physics and more recently also for the description of metallic clusters. As a detailed 
discussion is presented in the quoted papers we here focus our attention on the aspect of interest in the present 
context. 

The single particle potential 


Vig,z) 




2z^ - 3zg\ 


(3) 
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is, for 6 > 1 (6 < 1), a quadrupole deformed harmonic oscillator of prolate (oblate) shape with an additional octupole 
term; in fact the term multiplying A is proportional to r^P 3 (cos 0 ) with P 3 the third order Legendre polynomial. We 
use cylindrical coordinates 2: and q = ^■ For A 7 ^ 0 this is a two degrees of freedom system which is non- 
integrable. For 6 < 1 it turns out that the switching on of the octupole term very quickly gives rise to classically chaotic 
behaviour while the onset of chaos is barely discernable when b > 2. The statistical analyses of the correpsonding 
quantum spectra reveal the expected results in that for & < 1 the Wigner surmise is obtained for the nearest neighbour 
distribution while for 6 > 2 a Poisson distribution is obtained, and for particular values of A even a new shell structure 
emerges. The latter is understood in terms of corresponding classical periodic orbits [Qj^. 

In the spirit of the present paper the quantum mechanical findings should be directly obtained from the matrix 
structure of the associated Hamiltonian and the exceptional points related to it. The appropriate basis where = 
p^/2m + + z^lb‘^)l2 is diagonal is given by the occupation numbers n_L and riz. Note that the z-component 

of the angular momentum is conserved and we consider here only = 0. The arrangement of the quantum numbers 
in the matrix j is determined by the arithmetic ascending order of the unperturbed levels The 

arrangement will therefore depend on the value of b. The selection rules restrict entries in to n'^ = riz ± {2k + 1) 
and n'j_ = n±zt {2k) with fc = 0,1, 2,.... It turns out that for b>2 the entries vanish in the blocks of which refer 
to the subspaces in which the unperturbed energies are degenerate. As a consequence, the perturbative expansion 
at A = 0 of the eigenvalues starts with the quadratic term. From the discussion in the previous section it follows that 
a great number of exceptional points is trapped at A = 0 and we therefore expect the onset of quantum chaos only for 
values of A at an appreciable distance from zero. Using 544 x 544 matrices we illustrate in Fig. (4a) the distribution 
of the real parts of the exceptional points for the matrix problem + XH^ with 6=2. The bulk of the exceptional 
points occur in fact at values of A which fall outside the physical range (for A/Ac > I the potential no longer binds). 
In Fig.(4b) we diplay the corresponding distribution for 6 = 1/2. Now we obtain the maximum density of exceptional 
points as soon as the parameter A is switched on. This is in accordance with the results of the previous section since 
now the linear term does occur in the expansion of each level around A = 0 as there are non-zero entries in the blocks 
of which refer to the degenerate subspaces. The effect of this difference in the distribution of the exceptional 
points manifests itself in the different behaviour of the respective quantum spectra in that a Wigner distribution is 
found for A > 0 when 6=1/2 while for 6 = 2 a Wigner distribution never develops for A < Ac. 

V. SUMMARY AND DISCUSSION 

The previous section provides a fine example where the statistical properties of the quantum spectrum can be 
predicted from the properties of the individual matrices and alone. To obtain the distribution of the exceptional 
points as illustrated in Figs.(4) it is not necessary to solve the full problem -|- XH^ let alone to determine their 
positions (which would be impossible for more than 295,000 exceptional points). The distribution was found using 
the method developed in and explained in the Appendix. There, only some simple properties of and are 
used. While this result alone yields just another confirmation of the method employed in g], the new aspect of this 
paper lies in the prediction that, if the unperturbed levels are degenerate to second order, level statistics ascribed 
to quantum chaos are substantially suppressed initially and become manifest only for sufficiently large values of the 
perturbation. The argument comes about rather naturally from the behaviour of the exceptional points which in turn 
determine the degree of quantum chaos |^. While the previous section provides the physical relevance of our findings, 
their universal character is argued and demonstrated in section two and three. To the best of our knowledge, this is 
the first example where even finer details with regard to level statistics can be extracted from the distribution and 
the general behaviour of the exceptional points. We believe that more refinement can eventually even predict new 
shell structure as the one discovered in the model discussed. Work towards this aim is in progress. 


APPENDIX A: DISTRIBUTION OF THE REAL PARTS OF THE EXCEPTIONAL POINTS 

We use the method of the unperturbed curves [Q, where the actual spectrum is approximated by simple algebraic 
curves. The intesection points of the approximate curves are then identified with the real parts of the exceptional 
points. Since the aim is to obtain the distribution rather than the excat positions, the approximate curves suffice. 
They are obtained from the individual matrices and only. The requirement that the actual spectrum and the 
approximate curves coincide exactly for small and for large values of A leads for the prolate case (6 = 2) to the form 

( £i + 7 zA^ for A < 0.9Ac 

m) = ,_ 

l ^ujfX^ + SojfaiX^ + CiX + di for A > 0.9Ac 
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where the £i and uii are the eigenvalues of and , respectively, and 


£i- £n 

with U being the orthogonal matrix that diagonalises . The coefficients Ci and di are determined so as to smoothly 
match the two curves at A = 0.9Ac. We note that the second order perturbation is correct up to third order terms in 
A and is in fact very near to the actual spectrum for A < 0.9Ac, while the expression for A > 0.9Ac is correct up to 
terms of the order 1/A. 

For the oblate case {h = 1/2) the spectrum is, for the purpose considered, sufficiently well approximated by the 
unperturbed lines 


Gi{X) — £i-\- XcOi- 

Recall that the eigenvalues £i and Wj are different from those in the prolate case. 
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Figure captions 


Fig.l Exceptional points of the 9-dimensional model in the complex A-plane for /i = 0 (a) and /i = 1/2 (diamonds 
in (b)) and /r = 1 (crosses in (b)). Note the difference of the scale in (a) and (b) and the motion away from the origin 
when ^ increases. 

Fig.2 Spectra of the 9-dimensional model for ^ = 0 (a) and /i = 1 (b). For better illustration a partial view is 
given for ^ = 1/2 (c) and ^ = 1 (d). Note how the position of the avoided crossing not only moves outwards but also 
becomes more pronounced with increasing value of /r. 

Fig. 3 Sample averages of nearest neighbour distributions of 304 levels for /i = 0 (a) with an eightfold degeneracy 
at A = 0 and for /r = 1 (b); both distributions are for A = 0.2. 

Fig.4 Distribution of the real parts of the exceptional points for the prolate case (a) and the oblate case (b). 
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